Vortex studies in superconducting Ba(Feo.93Co .o7)2As 2 
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We present small-angle neutron scattering (SANS) and Bitter decoration studies of the supercon- 
ducting vortices in Ba(Feo.93Coo.o7)2As2. A highly disordered vortex configuration is observed at all 
measured fields, and is attributed to strong pinning. This conclusion is supported by the absence of 
a Meissner rim in decoration images obtained close to the sample edge. The field dependence of the 
magnitude of the SANS scattering vector indicates vortex lattice domains of (distorted) hexagonal 
symmetry, consistent with the decoration images which show primarily 6-fold coordinated vortex 
domains. An analysis of the scattered intensity shows that this decreases much more rapidly than 
expected from estimates of the upper critical field, consistent with the large degree of disorder. 

PACS numbers: 74.25.Qt, 74.70.Dd, 61.05.fg 



The recent discovery of superconductivity in 
LaFeAsOf 1 - with critical temperatures increasing to 
38 K upon doping or 43 K when subjected to hydro- 
static pressure,—^ has sparked a strong interst in these 
materials. However, synthesis, and especially the growth 
of single crystals, of these materials have proven to be 
a challenge. It was therefore a significant step forward 
when it was reported that the intermetallic compound 
BaFe2As2 could be rendered superconducting by doping 
either K on the Ba-site^^ or Co on the Fe-site^i 7 - This 
discovery enabled the possibility of obtaining large, 
high-quality single crystals grown from flux. 

Here we report on combined small-angle neutron 
scattering (SANS) and high resolution Bitter dec- 
oration studies of the vortices in superconducting 
Ba(Fe .93Coo.o7)2As2^ The SANS measurements show a 
ring of scattering, indicating a highly disordered vortex 
configuration, and with a rocking curve extending be- 
yond the measurable range. This result is confirmed and 
extended to lower fields by the Bitter decoration images. 
The disordering is attributed to strong vortex pinning, 
which is supported by the absence of a Meissner rim in 
decoration images obtained close to the sample edge. The 
field dependence of the magnitude of the SANS scatter- 
ing vector indicates vortex lattice domains of (distorted) 
hexagonal symmetry. An analysis of the scattered in- 
tensity due to the vortices shows a rapid decrease with 
increasing applied magnetic field, significantly exceeding 
what would be expected based on estimates of the upper 
critical field. 

Small-angle neutron scattering (SANS) measure- 
ments were performed on a 260 mg single crystal of 
Ba(Feo.93Coo.o7)2As2- The single crystal used in the 
SANS experiment was grown from a FeAs/CoAs mix- 
ture, and had a critical temperature T c = 21 K with a 
relatively narrow transition width AT C = 2 K measured 
by dc magnetization^ The experiment was performed at 
the D22 SANS instrument at the Institut Laue-Langevin. 



Incident neutrons with wavelength A„ = 6 — 14 A and 
wavelength spread of AA ra /A„ = 10% were used, and the 
diffraction pattern collected by a position sensitive detec- 
tor. Measurements were performed at 2 K in horizontal 
magnetic fields up to 810 mT, applied at angle of 5 — 10° 
to the c axis to reduce background scattering from crys- 
tallographic defects. All vortex measurements were done 
at 2 K, following a field cooling procedure. Background 
measurements obtained at 25 K > T c were subtracted 
from the data. 

The high resolution Bitter decoration experiments 
were performed on a Ba(Feo.93Coo.o7)2As2 single crystal, 
on an as grown surface perpendicular to the c axis, for 
applied magnetic fields between 1 and 32 mT along the 
c axis and following a field cooling procedure. The crys- 
tal surface was flat and shiny except in areas with growth 
steps. The initial sample temperature was 1.7 K, increas- 
ing to 5 - 6 K at the end of the decoration process. The 
sample was subsequently imaged using a scanning elec- 
tron microscope (SEM) in the second electron emission 
regime, to locate the islands of iron particles which had 
decorated the sample at the vortex positions^ 

Fig. QJa-d) shows vortex diffraction patterns obtained 
at four different applied magnetic fields between 810 and 
450 mT. For all the magnetic fields investigated, a dis- 
ordered vortex "lattice" (VL) is observed, evident by a 
"powder diffraction" ring of scattering instead of well de- 
fined Bragg peaks. It is reasonable to suppose that the 
disordering of the VL is caused by pinning of vortices to 
defects in the sample. Nonetheless, the ring of scatter- 
ing indicates that short range order persists, with small 
ordered VL domains randomly oriented with respect to 
one another. The intensity variation around the rings 
of scattering suggest the existence of maxima along the 
horizontal and vertical axis corresponding to the {100} 
directions, indicating a VL preferred orientation even if 
it insufficient to cause the coalignment of all domains. 

Scattering from the vortices is expected at a radial 
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FIG. 1: (Color online) Vortex imaging in Ba(Feo.93Coo.o7)2As2 by SANS (a-d) and Bitter decoration (e, f). The SANS vortex 
diffraction patterns were obtained at 2 K and applied fields of 810 (a), 740 (b), 630 (c) and 450 mT (d), following a field cooling 
procedure. The field was applied at a 5 — 10° angle with respect to the crystalline c axis to reduce background scattering 
from crystallographic defects. For all four fields data was collected at a single (centered) orientation of the magnetic field with 
respect to the neutron beam direction, and background measurements obtained at 25 K > T c were subtracted. The data is 
smoothed and the central part of the detector is masked off. The decoration images were obtained in a magnetic field of 4 mT 
in the center of the sample (e), and in the vicinity of the sample edge (f). Small ordered domains can be observed as indicated 
in (e). The inset to (e) show the Fourier transform of the entire decoration image. 



distance from the center of the detector which is propor- 
tional to the VL scattering vector, q. As the magnetic 
field, and hence the vortex density decreases, the average 
intervortex spacing increases as a cx (<fo/ B) 1 / 2 , where 
4> Q = h/2e = 20.7 x 10 4 TA is the flux quantum. Low- 
ering the field we therefore expect the ring of scattering 
to shrink, since q oc 2n/a tx (B/ifio) 1 / 2 decreases. This is 
clearly seen in Fig. [IJ a-d), confirming that the scattering 
is indeed due to the superconducting vortices. 

Results of Bitter decoration measurements obtained in 
a magnetic field 4 mT are shown in Fig. [TJc, f). A close 
inspection of the decoration image in Fig. [IJe) reveals 
primarily 6-fold coordinated vortices, as well as the ex- 
istence of small, well ordered domains with a size of ~ 5 
vortex spacings. However, no long-range orientational or- 
der is established, as evident from the Fourier transform 
of the decoration pattern image shown in the insert. The 
Fourier transform is directly comparable to the diffrac- 
tion patterns in Fig. [IJa-d). Combined with the SANS 
results this confirms that the disorder is predominantly 
orientational and longitudenal (along the vorticcs/ficld 



direction), but not Ad-disorder (vortex spacing). Simi- 
lar results were obtained at fields as low as 1 mT and as 
high as 32 mT. Computing the average area associated 
with each vortex in a 1 mT decoration image (not shown) 

we find a magnetic flux per vortex of (21 ± 1) x 10 4 TA , 
in agreement with 4>q. Further proof of strong bulk pin- 
ning is the absence of a Mcissner rim (low vortex density 
region) in the decoration image shown in Fig.[ljf), which 
was obtained near the edge of the sampled 

Our observation of a disordered vortex arrangement 
is consistent with recent reports of STS imaging in 
Ba(Feo.93Coo.o7)2As2^ However, compared to small- 
scale direct space imaging like STS, our results demon- 
strate that the disordering is a bulk effect, not restricted 
to the sample surface or isolated regions. 

The field dependence of the scattering vector is also ev- 
ident from Fig. [21 which show the radial intensity distri- 
bution corresponding to each of the diffraction patterns 
in Fig. QJa-d). Fitting each radial intensity curve to a 
Gaussian it is possible to obtain a quantatitive measure 
of both the magnitude of the scattering vector and the 
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FIG. 2: (Color online) Radial intensity distribution for the 
four diffraction patterns in Fig.[TJa-d). The curves have been 
offset by 2 Cts./Std. Mon. for clarity. Also shown is a Gaus- 
sian fit to the data at each field. 



scattered intensity. Furthermore the FWHM of the fits 
yield Aq/q ss 0.3, which is ~ 50% larger than the exper- 
imental resolution and consistent with scattering from 
small ordered domains of size ~ (irAq/q) -1 , correspond- 
ing to a few VL spacings in agreement with the decora- 
tion results. If one considers a rhombic VL unit cell as 
shown in the insets to Fig. [3^b), the scattering vector is 
given by q — 2n(B/cj)Q sin/3) 1 / 2 where (3 is the opening 
angle. Here (3 = 60° (or 120°) corresponds to a hexago- 
nal VL and f3 = 90° to a square VL. Fig. Ela) shows the 
measured scattering vector as a function of applied mag- 
netic field, compared to that expected for respectively 
a square and a hexagonal VL. With the exception of a 
single data point, all the measured scattering vectors fall 
between (7h cx and q sq , although it is closer to the former. 
This is consistent with the observation of primarily 6-fold 
coordinated vortices in the decoration images. Fig. [3Jb) 
shows the measured scattering vector normalized to q sq . 
Within the error bars, no field dependence of the normal- 
ized scattering vector is observed over the measured field 
range. 

While a hexagonal VL is expected for an isotropic ma- 
terial, it is also well known that most tetragonal su- 
perconductors posses enough anisotropy to stabilise a 
square VL above a certain magnetic field i 11 ' 12 ^ 3 ' 14 ^ 5 
The origin of the anisotropy can arise, either from the 
Fermi surface of a given material coupled with non-local 
electrodynamics^ or it can be due to anisotropic pairing 
in the superconducting stated Given that a near hexag- 
onal VL symmetry is observed in Ba(Fco. 93000.07)2^82 
at all the measured fields, it is tempting to conclude that 
the supercondcuting state in the basal plane of this mate- 
rial is highly isotropic. However, it has been shown that, 
at least in the case of a square VL stabilised by a Fermi 
surface anisotropy, the critical field for the transition to 
the square symmetry depends strongly on the so-called 
non-locality range, which decreases with increasing levels 
of impurities (decreasing mean free path) ^ From our re- 
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FIG. 3: (Color online) Comparison of the measured scatter- 
ing vector (a), and the ratio q/q aq (b), to that expected for 
respectively a square (dashed lines) and a hexagonal (dot- 
dashed lines) vortex lattice. The solid line in panel (b) is the 
average g/q sq = 1.055. The different symbols corresponds to 
measurements performed at different neutron wavelength, X n . 



suits it is not possible to conclude which mechanism is re- 
sponsible for the hexagonal VL in Ba(Feo.93Coo.o7)2As2. 

We now turn to the integrated scattered intensity, ob- 
tained from Gaussian fits to the radial intensity distri- 
bution as shown in Fig. [2] From a measurement of the 
absolute reflectivity, R, it is possible to determine the VL 
form factor by 



27T7 2 A 2 t 



R, 



(1) 



where 7 = 1.91 is the neutron gyromagnetic ratio and 
t is the sample thickness However, to put the form 
factor on an absolute scale, the scattered intensity from 
the vortices should be integrated as the sample and mag- 
netic field is rotated with respect to the incoming neutron 
beam, causing the reflections to cut through the Ewald 
sphere. In the case of Ba(Feo.93Coo.o7)2As2 scattering 
from the vortices showed no measurable change in inten- 
sity even for rotations of several degrees. This indicates 
that the high degree of disorder found in the basal plane 
is also present along the direction of the vortices, causing 
them to deviate substantially from the applied field direc- 
tion. In this regard the results in Ba(Feo.93Coo.o7)2As2 
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FIG. 4: (Color online) Field dependence of the vortex scat- 
tered intensity. The different symbols corresponds to mea- 
surements performed at different neutron wavelength, A n . 
The solid line is an exponential fit to the data with a slope 
of —0.82 (1/T). The dashed and dot-dashed lines show the 
form factor calculated using the London model (Eq. (J3}) with 
c = 2 and and 1/2, and H c2 = 40 T. 



are reminiscent of those obtained on G&Cq^ Our results 
are consistent with the report of Prozorov et al. of a 
highly disordered vortex phase in Ba(Feo.93Coo.o7)2As2, 
based on the analysis of transport measurements 3i 

Our inability to measure a proper rocking curve means 
that a determination of the form factor on an absolute 
scale is not possible. On the other hand, if we assume 
a constant rocking curve width, the radial integrated in- 
tensity will be proportional to |F(<7)| 2 . Due to the strong 
pinning in the sample, and relatively small field range 
investigated, we we believe that this is a resonable as- 
sumption. The form factors obtained in this fashion are 
shown in Fig. 0] on a logarithmic scale. Here it is impor- 
tant to stress that even though the absolute magnitude 
of |.F(g)| 2 is not known, the relative field dependence can 
still be determined. The data are well fitted by an ex- 
ponential field dependence as shown by the solid line in 
Fig.H 

Several models exist for the form factor field depen- 
dence. By far the simplest is based on the London model, 
extended by a Gaussian cut-off to take into account the 
finite extent of the vortex cores: 



F(q) 



B 



(2) 



1 + (Xq) 2 

Here A and £ are respectively the penetration depth and 



coherence length, and the constant c is typically taken to 
be between 1 /4 and 2.— Even though more realistic mod- 
els for the form factor are available, the London model 
is deemed adequate for the following discussion since the 
SANS measurements were performed at T ~ 0.1 T c and 
H <C H C 2- For simplicity we will also assume a square 
vortex lattice for the following analysis. It should be 
stressed that neither of these choices affect the following 
discussion in any significant fashion. For all fields studied 
(qX) 2 ;>> 1 and since q 2 = (2ir) 2 B/cj)o the form factor in 
Eq. ^ decreases exponentially with increasing magnetic 
field. 

As shown in Fig. |4j the measured intensity is indeed 
well fitted by an exponential decrease. If we use the 
fitted slope to estimate the superconducting coherence 
length we obtain £ = 10 nm. Here we have used c = 0.5 
which in the past been found to yield reasonable val- 
ues for £ in other materials^ However, this value of £ 
is roughly 3.5 times higher than expected from the ex- 
trapolated upper critical field of ~ 40 T at 2 Kj2 For 
comparison the H C 2 corresponding to this "measured" 
£ is only 3.3 T, more than an order of magnitude less 
than the measured upper critical field. The difference is 
emphasized in Fig. [4] which shows the range of the cal- 
culated form factor (shaded area) corresponding to the 
real H C 2 and the two extreme values of c. From this it is 
clear that the measured form factor falls off much faster 
than expected. We believe that this rapid decrease in 
the scattered intensity is due to the high degree of disor- 
der of the "vortex lattice" in Ba(Feo. 93000.07)2 As2- This 
result is indicative of the existence of a vortex glass or 
randomly oriented domains of a Bragg glass, consistent 
with the results of Klein et al£^ 

In summary, we have performed SANS and Bitter 
decoration studies of the vortices in superconducting 
Ba(Feo.93Coo.o7)2As2. At all the fields investigated, the 
results show a highly disordered vortex configuration was 
observed, most likely due to strong pinning in this mate- 
rial. Further measurements at higher fields and temper- 
atures should be carried out, to fully explore the vortex 
phase diagram in this and related comppounds. 
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